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Abstract. We formulate the problem of the Fermi Edge Singularity in non- 
equilibrium states of a Fermi g matrix Riemann-Hilbert problem with an 
integrable kernel. This formulation is the most siiitable for studying the singular 
behavior at each edge of non-equilibrium Fermi states by means of the method 
of steepest descent, and also reveals the integrable structure of the problem. We 
supplement this result by extending the familiar approach to the problem of the Fermi 
Edge Singularity via the bosonic representation of the electronic operators to non- 
equilibrium settings. It provides a compact way to extract the leading asymptotes. 
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1. Introduction 

The FES (Fermi edge singularity) [H [21 H] is observed in absorption of X-rays in 
metals as a power law peak at the Fermi Edge of a degenerate Fermi gas. There a 
sudden removal of a localized electron from a hard core atomic shell creates a potential 
which disturbs the electronic gas, thus producing a power low spectrum of electronic soft 
modes. In recent years the FES has also been demonstrated in tunneling experiments 
[3 El m E]. There a single electron can change the capacity of a contact producing a 
similar disturbance to the electronic g clS clS db localized hole. As a result, one observes 
a power law in tunneling current vs. the bias voltage: I{V) ~ y-'ia+ka'^ where 
6 = ira is the scattering phase of the ensuing potential and k is the number of scattering 
channels. In the case of an attractive potential (a > 0) the current peaks at the Fermi 
edge. 

One of the reasons of interest in the FES, and our own motivation in studying it, 
is that the origin of the FES may be found to be ascribable solely to Fermi statistics. It 
has been studied over at least five decades, has been well understood, and is considered 
as one of the fundamental quantum phenomena in electronic physics. Early theoretical 
papers [H El [3] on the FES and the related phenomenon of Orthogonality Catastrophe 
[To] were proved to be influential well beyond FES. They are at the foundation of the 
modern physics of electronic systems in low dimensions. 

In a degenerate Fermi gas anything but the leading power asymptotes is rarely of 
any interest. The reason for that is that everything else, except the leading power, 
depends on details on band structure, tunneling contacts, etc., and lacks of universal 
character. 

A different situation occurs in a non-equilibrium Fermi gas. There the energy scale 
of non-equilibrium features can be much smaller than the Fermi scale and can be seen 
in the spectrum of absorption or tunneling. In this case, additional features become 
universal. 
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Figure 1. Occupation of electronic states in structured Fermi sea. 

As a prototype of a non-equilibrium state we consider a state where the Fermi 
distribution consists of steps at Ei < E2 < . . . -^2^+1 such that no states are occupied 
in the interval between energies E2i-i,E2i, where i = 1, . . . ,n. We denote by Eq the 
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bottom of the filled Fermi sea, which is assumed to be far away Fig. [T] Assuming that 
the scattering phase does not change within a wide range of the conducting band, say, 
between £"0 = to = A the spectrum will be a transcendental universal function of 
V/Eij, where Eij = Ei — Ej. 

These structured non-equilibrium states described above inevitably appear in the 
evolution of an arbitrary semiclassical Fermi state [11] . They were also realized in some 
nanoscale devices (see e.g. [12j). 

In a non-equilibrium state the absorption spectrum (or tunneling current) is 
a transcendental function and elementary methods can determine its asymptotic 
expansion in various regimes. In this paper we develop a framework aiming to 
characterize the universal part of the spectrum. We derive a 2 x 2 matrix Riemann- 
Hilbert problem derived from an inversion problem of an integrable kernel. A similar 
matrix RH-problems appeared in studies of various fermionic correlation functions 
and correlation functions of eigenvalues of Random Matrices [13]. One of a proven 
advantage of this formulation is that RH-problem is the most suitable for application 
of the steepest descent method. The latter prompts the leading asymptotes at Fermi 
edges. Besides the fact that the RH-problem connects FES to a number of physically 
unrelated but mathematically equivalent (often well studied) problems, we think that 
these links and analytical apparatus they bring to physics of the FES are important for 
a deeper understanding quantum non-equilibrium phenomena. A somewhat alternative, 
but essentially equivalent approach is to establish a set of non-linear integrable equations 
[14j . We do not discuss this approach here. 

As far as the main asymptote is concerned, we will show that in the setting described 
above the tunneling current at voltage close to a Fermi edge {eV — Ei) <^ |-Ej,j±i| 
reads: 

/(V) oc Ai|ey-Ei|('"-')"'-'"'", (1) 

2n 2n 

A, = const n ^z'^"^'""' n Ei^^'"^ (2) 




Figure 2. A schematic plot of tunneling current for small a < (left panel) and 
a > (right panel), solid lines show computed power law asymptotes. Dashed lines 
interpolate between resonances. 
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where Ei = ±1, if Ei is an lower/upper edge of an occupied band, respectively. 

Some noticeable features of this result are: (i) the exponent depends on the total 
number of bands, but stays the same for all upper (lower) edges of occupied intervals 

= (2na^ — (a — — Ei/e\~^; (ii) if the potential is attractive a > (which is 

the common case) the current features a peak with a power low decay at upper edges 
towards increasing energy; and in contrast the current is suppressed at lower edges; if 
the potential is repulsive, a < 0, and sufficiently small, a peak appears at the edge E2k 
with a power law to the left to the edge Fig|2] (iii) the logarithm of the factor Ai can be 
treated as the electrostatic energy of particles with alternating charges ±a positioned 
at the upper/lower edges with an insertion of a unit charge at the edge Ei. 

Fermi Edge singularities with a structured Fermi distribution have been studied 
before. In 1984 Combescot and Tanguy [13 [ISl [IZ] considered a situation where 
the interval {Ei,E2) of the band (0,A) is occupied, while the interval {0,Ei) starting 
from the bottom of the band is unoccupied. Later 2004 the Abanin and Levitov [18] 
considered FES with a two steps Fermi distribution. A more general situation has been 
considered in the recent papers by Gutman et al [19] , and by the authors [11] , where Eq. 
([T]) have been obtained. The approach employed in these papers (with the exception 
of Ref.[lT]) is based on the expression of the tunneling current as a determinant of a 
Fredholm operator. A basis of particle-hole excitations of the structured Fermi sea has 
been used to write the Fredholm operator in early papers [151 HBl [E] ■ This approach 
has been developed by Othaka and Tanebe in 1984. Contrary one-particle basis 
in an empty vacuum has been used in later articles [iHl [19]. Naturally the basis of of 
particle- hole excitations of the Fermi sea, employed in earlier papers [ISl [161 [T71 H] , which 
already captures the many-body physics of the Fermi step forward in obtaining 

the desired results. We employ this approach here. In addition, for illustrative purposes 
(Appendix Appendix B ), we also present a simple a compact method to capture a leading 
singularity developed in Ref.[lT]. That method is based on the bosonic representation 
of electronic operators in a Fermi state with multiple edges. 



2. Tunneling current. 

The tunneling setting of FES is as follows: a Fermi gas is in contact with a localized 
resonant level (a quantum dot). It is initially uncharged and provides no scattering to 
electrons. When an electron tunnels to the dot, it suddenly charges the dot, switching- 
on a small potential H ^ H' = H + U localized at the dot [S] . We assume no further 
interaction, no dissipation, and we ignore spin and channels. 

The tunneling current is given by the golden rule [H [H] • In units of the tunneling 
amplitude, I{uj)\hu}=ev+Ei reads 

I{uj) oc Re r e''^^''-'^^G{ti,t2)dT, (3) 
Jo 

G{ti, t2) = {n\e'^''ce'^'^''-''^c^ e-'^''\n) (4) 

Here c = X^eQ is an electronic operator in the position of the dot, assumed to be the 
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origin, while is an electronic mode with energy e and \Q) is the structured Fermi state 
of interest. We also assume that the bias voltage is with respect to Ei. 

Next we assume that the potential is regular within a wide energy range which 
exceeds the energy range of features of the structured Fermi state and that the tunneling 
time is sufficiently small. Under this assumption the energy dependence of the scattering 
phase 6 caused by potential U can be dropped. This amounts to a downwards shift of 
the energy levels by a constant amount a (in units of level spacing): e — ?► e — a. 

In Ref. [3] it has been shown that the vertex operator e"*^ implements a shift 
of momenta: the perturbed Hamiltonian and perturbed states are seen as similarity 
transformation of the unperturbed ones H' = e~°-''^He"''^ and = e°''^\Q). Here an 
operator cp is a chiral canonical Bose field related to the chiral part of electronic density 

fit) = ^ E e^'Ve/e, P. = E 4ce+e- (5) 

Then Green's function reads 

G{t,,t,) = {n\c{t2)e-'^^^''^e'^^^'-^^c^{t^)\Q), (6) 
where c(t) = J2e es'^*^. This formula is standard. 



3. Fredholm Determinants 



Following [T], Green's function can be understood as consisting of three multiplicative 
factors G{ti,t2) = \ {Q'\Q)\'^e'" ■ L - an overall normalization: 

where = e^'^^^^lQ), (7) 



closed loops 

c _ 

" m'\n)\ 



^ - i/n/in\i2 ' (8) 



and open lines 

The two latter objects can be cast in the form of Fredholm determinant by means of 
the Wick theorem. We remind the major formulas. 

Consider coherent states of Gl{oo). These states are obtained by transforming the 
ground state of the Fermi gas \ A) = g{A)\0) by an exponent of a bilinear form of Fermi 
operators g{A) = e^";'^'''''^"'^'', where A^^ri is an arbitrary gl{oo) matrix. Our structured 
Fermi state is a coherent state. 

For arbitrary coherent states Ai, A2 and arbitrary (7/(00) matrices B, G the following 
holds 

{AMB)giG)\A2) 
{AMB)\A2){AMC)\A2) 

where 



det(l + K), (10) 



K,^,,,=J2M,^,MuB,As)MlJAs,G,A2) (11) 
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M,^^{A^,B,A,) = {AMB)clc^\A^) (12) 

are matrix elements of operator g{B) between the state (^^(v4)| and a state cj^c^l^a), 
where a particle-hole pair is added to an arbitrary chosen coherent state l^s). The 
result does not depend of the choice of ^3. 

In order to obtain this formula one inserts a superposition of an arbitrary number of 
particle-hole excitations into a chosen coherent state |y43)(y43|, apply the Wick theorem 
to each term and sum them up. 

The variation of (10) 51ogdet(l + K) = tr [(1 + K)^-'^5K] gives another known 
formula 

{A,\clMB)9{C)c.M2) 



{A\g{B)g{C)\A') 

where 



tr 



l + K)-ip , (13) 



P.,,(ei,e2) = M,,r,{A,,B,A2)M;^JA,,C,A2), (14) 

and M^rii^i, B, A2) = {Ai\ceg{B)clj\A2) . Specification of these formulas: (y4i| = {A2\ = 
{n], g{B) = g^{C) = e"'^'^'-*-' prompts a determinantal representation of the current 

e^ = det(l + K), L = tr[(l + K)-^P] (15) 

with a kernel 

7^(61,62) =E^e„,(tl)M;,^(t2), (16) 
V 

P,,, = M,{t^)M;{t2), (17) 

where 

_ {n\e-^-(Mcv\^) .... 
M.,,{t) , (18) 

^'^^^ WW) — ^^^^ 

are matrix elements of the vertex operator between states where a particle-hole pair, or 
just one particle are added to the state [|] These are general formulas valid for any 
coherent state \ fl). 

The formulas is further specified since the structured Fermi state we are considering 
is an eigenvalue of the Hamiltonian. In this case 

M^ri is an overlap between the state which appears after the shake-up and a particle- 
hole excitation of the state \Q) before shake-up with energy e ^ Q and rj & Q, where we 



I Formulas equivalent to ( 16 1 for the ground state can be found in [3] . 
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denote = U"ro'^(i?2i, -^24+1) as the set of occupied single particle states in {Vl,e) 
and [Q] e, rj) are states where an extra particle or a particle-hole pair is added into 1^2). 
Similarly is the overlap of states \Q) and the state with an added particle with 
energy e. 



4. Matrix elements 



The following formula helps evaluating the matrix elements ( 20pT ): If |e') is a single 



particle eigenstate of the perturbed Hamiltonian and |e) is a one-particle eigenstate of 
an unperturbed state then their overlap (in units of level spacing) is (e'|e) = f^r^- An 
extension of this formula to many particle states (e| = (ei, . . . , eAr| and |e') = \e[, . . . , e'^) 
gives a Cauchy determinant 



^ \^ lAjy 1 ni>j(ei-ei)(ei-4> 



= det = ' ^ \ " . (22) 



sin(7ra)y ^i^^'j riijlq - 



With the help of this formula, the matrix elements in ( 16 ) can be computed in a manner 
similar to Ref. [1]. While computing one must take into account that the set of occupied 
single particle levels in Vt and Vt' are shifted by a with respect to each other. Then the 
problem is reduced to an electrostatic problem of placing a dipole or a charge into a 
Coulomb plasma confined in the intervals U"^]^(i?2j, E2i+i)- It gives the overlap between 
states {VL\ and \VL') generalizing Orthogonality Catastrophe formula pD]. Up to an Ei 
-independent constant factor it reads 

i>j 

The results for matrix elements are 

^ g|(,-e)i^le)fM, = e^V(e), (24) 
e — rj 

2n 

r(6) = n (e - i?.)''\ ein- (25) 

2n 

s{r^) = '—X{{E,-r^r^\ 

1=1 

where e ^ Q, r] E Q are energies of particles and holes. A short sketch of these 



calculations is found in Appendix Appendix A 
Summing up, the kernel reads 

K{e„e,) = et(-*--*^)r(6i)r(e2)^^^^^^^, (26) 

^1 — ^2 



,i.r^M^^ (27) 

5. Integrable kernel 



The next step is to invert the Fredholm kernel K. It can be done in a straightforward 
manner similar to [U [T51 HSl HZ] employing the Wiener-Hopf method at every edge. 
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However, calculations become more structured if we use the integrable property of the 
kernel. Integrability is general property of free fermion correlators (see [201 

A kernel is called integrable if it has the form 

^1 ~ ^2 a=l 



In the case of a structured Fermi sea / = 2, and, as follows from ( T6p4 ): 



g^ie) = g(e, r)^72(e), /2(e) = -Qie, r)/i(e), (28) 
g,(e) = e^*V(e), A(e) = e-*^*V(6), 

Let F = (1 + K)~V (we denote / = (/i,/2)) be a solution of the singular integral 
equation 

F(ei) + / ir(ei, 62)^(62)^62 = /(6i), 61 ^ n. (29) 



The time derivative of closed loops contribution and a contribution of open lines (15) 
are expressed through the solutions F = (Fi, F2) 

f ^t..((l.K)-.f)4/^^tef.-,.F.)* (30) 

g2Fide, (31) 

6. Matrix Riemann-Hilbert problem 



The Fredholm equation (29) is sufficient to obtain the singular behavior at Fermi edges. 
However, it is instructive to cast the FES problem as a matrix RH problem along the 
lines described in [T3]. In that form, the FES problem falls in the general scheme of 
integrable problems. In addition, the RH-problem is the most suitable for analysis near 
edges [SU]. 

The central object of the RH problem is a matrix-valued functions m{e) analytic in 
a complex 6-plane cut along the unoccupied intervals U(i?2j-i, E2i) Fig. [3| defined such 
that at infinity m approaches the unit matrix, and that its boundary value on the cuts 
m± = m{e ± iO) connects vector the F to the vector / as 

F{e) = m^{e)f{e), e e Q. (32) 

In Ref. [13] it has been shown that the matrix is a solution of the RH-problem: 

m+v = m_, Vafi = Sai3 - 2mfagi3- (33) 
In the case of FES 

v{e) = l + 2me-'n^W\e)(^^, _M . (34) 



Eqs (30 34) constitute the matrix RH-problem for FES. As typical for other integrable 



RH-problems, a similarity transformation can be found to reduce the jump matrix v{e) 
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to a constant matrix, such that analytic behavior in the energy dependence of the kernel 
will be translated to into the analytic nature of singularities of the solution at infinity. 
We do not do this here. 

Being specified for a one-edge problem (n = 0), in units of upper and lower cut-offs, 

read 

r(o)(e) = e^'=^ ^^e > 0; (35) 

Q(°)(e,r)=f^)7 et-iz!^!!!!^^, (36) 
V vr / Jekv<o e — 1] 

where we count energy from the edge. In this case the RH-problem is solved by 

elementary means. In fact technically it easier to proceed directly through the integral 

equation following 

7. Method of steepest descent and the leading singularity 

The asymptotic behavior at the edges can be found by the steepest-descent method 
described in [20] . 



Figure 3. Unoccupied electronic states are between edges i?2fc-i and E2k, k = 1, . . .n. 
The matrix m of the RH-problem (331 jumps on segments of real axis corresponding 
to unoccupied states (solid line). The steepest descent contour goes down vertically 
(dashed line) in the lower half plane. 



The steepest-descent contour starts from Fermi edges and extends to ioo in the 
lower half-plane as is in Fig. [3j Along this contour the rapidly falling exponential factor 



e fi'^'^ in (34) suppresses the jump of the matrix m except at small segments near the 
Fermi edges, where r(e) and Q{e) are singular. At energy close to an edge, say, Ei we 
estimate r{Ei + e) ^ rW(e) and Q{E, + e, r) ^ et^'" ^j^'''"Q^°He), and 

the problem reduces by a similarity transformation to the one-edge problem 

m + e) = z^f'\e); m + e) = 2rV°He); 

m{Ei + e) = Zim^^\e)Zr^- 

where script (0) indicates entities of the one-edge problem, t = ti+t2, and 1 is the 2x2 
unit matrix. 

The reduction of the multi-edge problem to the one edge problem using the RH 
steepest descent method is what allows us to solve the problem. Indeed, the contribution 
of each edge is well known. In units of spacing and up to a constant factor they are 

Cf~r-"', Lf-e.r^'^-i. (38) 
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We note that Ohtaka and Tanebe [1] showed how (38) originally obtained in Refs. [HE] 
by different means) follow from the integral equation (29). 

The leading asymptote for the structured Fermi sea can be obtained by a 
combination of the one edge problem and a similarity transformation. The similarity 
transformation does not affect the contribution of closed loops. Therefore each edge 

„ „(0) 

contributes equally. Summing them up, we obtain e ~ * . 

The similarity transformation and (31) give the contribution of open lines as a sum 
of one edge open lines each weighted by its own amplitude 

2n+l 

L= Ee-^^'^n^?"%^°^ (39) 

i=l j^i 

Combining the normalization the closed loop and the open line we obtain: 

\i<j / m ni^m 

where a multiplicative constant depending only on cutoffs has been omitted. Fourier 
transforming, this formula prompts the main result 
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Appendix A. Compnutation of (23 24) 



We sketch the computations of matrix element (23). They are not much different from 
similar calculations for a single edge problem [1] . We start from a general formula (22) 
describing an overlap of arbitrary electronic states consisting of a finite number A of 
occupied levels. It can be seen as the exponent of the electrostatic interaction energy of 
log-interacting particles with charges of size -|-a at e Q's and —a at G — a. Thus 
we have a regular alternating pattern of charges, the negative charges at an offset of —a 
with respect to the positive ones, the patten starts at Eq and persists up to E2n+i, with 
gaps at [E2i-i, E2i], i = 1, . . . ,n. The electrostatic energy is: 

E log f ] -N\oga. (Al) 

This sum can be evaluated in terms of Barnes functions, but we are interested only in the 
limit of large number of levels. There we may think in terms of a density of dipoles with 
a polarization vector a. The energy of this system is equivalent to energy of charges eta 
placed at edge Ei. The latter (up to a-dependent constant) is Z]i<j ^j^iC^^ log \Ei — Ej\. 
The exponent of the energy yields ( [23| . 

The Configuration of charges corresponding to is a set of dipoles discussed 
above minus a charge +1 at e G f2 plus a charge —1 at rj ^ Q. The energy acquires the 
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interaction energy between charges at rj and e and charges ±a sitting at edges Ei. This 
addition is J2i ^iO^iJog \ e — Ei\— log \ ri — Ei\) plus an a dependent constant. The origin of 
the constant is the energy of dipoles in the vicinity of a hole at e. It is J^eten log • 
Assuming that rj is far from edges the contribution goes from levels close to e. It gives 



log ('^). All together it yields (24). Computation of M, is similar. In that case only 
interactions between the extra charge at e ^ f2 and residual charges at edges contribute. 
It does not incur a constant factor. 

Appendix B. Computation through bosonic representation 

The leading singularity can be understood using a bosonic formalism. First we separate 
fast oscillatory modes at each edge 

c{t,Xo) =J2e^^^%{t). 

i 

Then we represent slow modes through components of the Bose field dx^^i = itpi(t)ipi{t) 
as 

V'.oc(5,nE^y/v^'^\ (B.i) 

The Bose field Q is a sum of its components ip = J2k Vt- Components of the Bose 
field represent particle-holes excitations close at each edge. At Ti/t ^ Eij they 
can be treated as independent canonical Bose fields. Their variances Di{ti,t2) = 
— {(pi{t2) — ipi{ti))'^ \Q) are not difficult to compute. As follows form (5), D2k-i 
are sums of (coser — l)/e over all possible energy of a particle-hole excitations provided 
that a particle is placed in the "gap" (£'2^-1, E2k) ■ Similarly D2k is the sum over energies 
of a hole-particle excitations provided that a hole is placed to the band (£'2^, -E'2fc+i)- 
Computing these integrals at r ^ h/\Eij\ one obtains 

A,(r) = - log r + ^ e, log \Ei,\ (B.2) 



The time independent term in (B.2) explains the prefactor in (B.I): the correlator 



{'4'i{ti)ipi{t2)) oc ^ then also has to be obtained for the Bose field. This yields 



In the Bose representation, Green's function ([6]) is a sum of edge components 

^(^1,^2) = E(^*n^r''>^''''G.(ti,t2), (B.3) 



Computing this, we obtain Green's function (40). 
References 



[1] P. Nozieres and C. T. de Dominicis. Singularities in the X-Ray Absorption and Emission of Metals. 
III. One-Body Theory Exact Solution. Physical Review, 178:1097-1107, 1969. 



Fermi Edge Resonances in N on- equilibrium States of Fermi Gases 



12 



G. D. Mahan. Excitons in Metals: Infinite Hole Mass. Physical Review, 163:612-617, 1967. 

K. D. Schotte and U. Schotte. Tomonaga's Model and the Threshold Singularity of X-Ray Spectra 

of Metals. Physical Review, 182:479-482, 1969. 
K. Ohtaka and Y. Tanabe. Theory of the soft-x-ray edge problem in simple metals: historical 

survey and recent developments. Reviews of Modern Physics, 62:929-992, October 1990. 
A. K. Geim, P. G. Main, N. La Scala, Jr., L. Eaves, T. J. Foster, P. H. Beton, J. W. Sakai, F. W. 

Sheard, M. Henini, G. Hill, and M. A. Pate. Fermi-edge singularity in resonant tunneling. 

Physical Review Letters, 72:2061-2064, March 1994. 

D. H. Cobdcn and B. A. Miizykantskii. Finite- Temperature Fermi-Edge Singularity in Tunneling 
Studied Using Random Telegraph Signals. Physical Review Letters, 75:4274-4277, December 
1995. 

I. Hapkc-Wurst, U. Zeitler, H. Frahm, A. G. M. Janscn, R. J. Haug, and K. Pierz. Magnetic-field- 
induced singularities in spin-dependent tunneling through InAs quantum dots. Phys. Rev. B, 
62:12621-12624, November 2000. 

Y. N. Khanin, E. E. Vdovin, L. Eaves, I. A. Larkin, A. Patane, O. N. Makarovskii, and M. Henini. 
Magnetic-field-induced Fermi-edge singularity in the tunneling current through an InAs self- 
assembled quantum dot. Soviet Journal of Experimental and Theoretical Physics, 105:152-154, 
July 2007. 

K. A. Matvcev and A. I. Larkin. Interaction-induced threshold singularities in tunneling via 

localized levels. Phys. Rev. B, 46, 1992. 
P. W. Anderson. Infrared Gatastrophe in Fermi Gases with Local Scattering Potentials. Physical 

Review Letters, 18:1049-1051, June 1967. 

E. Bettelheim, Y. Kaplan, and P. B. Wiegmann. Gradient Gatastrophe and Fermi Edge Resonances 
in Fermi Gas. ArXiv e-prints/ 1011.1993, November 2010. 

S. de Pranceschi, R. Hanson, W. G. van der Wiel, J. M. Elzerman, J. J. Wijpkema, T. Fujisawa, 

S. Tarucha, and L. P. Kouwenhoven. Out-of-Equilibrium Kondo Effect in a Mesoscopic Device. 

Physical Review Letters, 89(15):156801-+, September 2002. 
P. Deift, A. Its, and X. Zhou. A Riemann-Hilbert approach to asymptotic problems arising in 

the theory of random matrix models, and also in the theory of integrable statistical mechanics. 

Ann. of Math., 146:149 235, 1997. 
E. Bettelheim, A. G. Abanov, and P. B. Wiegmann. FAST TRAGK GOMMUNICATION: 

Quantum hydrodynamics and nonlinear differential equations for degenerate Fermi gas. Journal 

of Physics A Mathematical General, 41:2003— October 2008. 
M. Combcscot and C. Tanguy. Absorption-edge singularities for a nonequilibrium Fermi sea. I. 

Second-order perturbation theory. Phys. Rev. B, 50:11484-11498, October 1994. 
G. Tanguy and M. Gombescot. Absorption-edge singularities for a nonequilibrium Fermi sea. II. 

Second-order diagrammatic expansion. Phys. Rev. B, 50:11499-11507, October 1994. 
G. Tanguy and M. Gombescot. Absorption-edge singularities for a nonequilibrium Fermi sea. III. 

Determinantal nonperturbative theory. Phys. Rev. B, 52:11698-11710, October 1995. 
D. A. Abanin and L. S. Levitov. Timablc Fcrmi-Edgc Resonance in an Open Quantum Dot. 

Physical Review Letters, 93(12):126802-+, September 2004. 
D. B. Gutman, Y. Gefen, and A. D. Mirlin. Non-equilibrium ID many-body problems and 

asymptotic properties of Toeplitz determinants. ArXiv e-prints, October 2010. 
P. Deift and X. Zhou. A steepest descent method for oscillatory Riemann-Hilbert problems. ArXiv 

Mathematics e-prints, December 1992. 
M. Jimbo, T. Miwa, Y. Mori, and M. Sato. Density matrix of an impenetrable Bose gas and the 

fifth Painlcve transcendent. Physica D Nonlinear Phenomena, 1:80-158, April 1980. 
A. R. Its, A. G. Izcrgin, V. E. Korepin, and N. A. Slavnov. Differential Equations for Quantum 

Gorrelation Functions. International Journal of Modern Physics B, 4:1003-1037, 1990. 



